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1.1  INTRODUCTION 

The  Ocean  Measurements  Program  (OMP)  of  the  Naval 
Ocean  Research  and  Development  Activity  (NORDA)  has  iden¬ 
tified  the  measurement  and  characterization  of  upper-ocean 
vertical  shear  as  a  high  priority  component  of  its  mission. 
In  support  of  this  objective,  SAI  has  recently  completed 
reviews  of  relevant  results  of  current  research  on  upper- 
ocean  vertical  shear.  Lambert  and  Patterson  (1980)  reviewed 
the  measurements  up  through  1979,  emphasizing  vertical 
structure  of  currents  and  shear.  Their  summary  included  a 
discussion  of  the  particular  relevance  of  shear  to  Navy 
problems.  Grabowski  and  Hebenstreit  (1981)  completed  a  com¬ 
prehensive  review  of  upper-ocean  environmental  models,  with 
heavy  emphasis  on  internal/inertial  wave  models. 

Following  these  review  efforts,  several  studies  of 
shear  were  initiated.  An  internal-wave  spectral  model  was 
used  to  calculate  shear  statistics,  and  an  error  analysis  was 
carried  out  to  estimate  the  utility  of  WKB  scaling  for 
extending  such  models  into  the  seasonal  thermocline 
(Grabowski,  1980).  A  data  analysis  effort  focusing  on 
vertical  profiles  of  shear,  presented  some  generalizations  of 
the  descriptions  of  the  shear  field  for  different  background 
fields  of  stratification  and  mean  currents  (Patterson,  et_ 
al . ,  1981).  During  the  past  year  SAI  has  carried  out  two 
parallel  efforts  to  describe  upper-ocean  shear.  A  modeling 
effort  has  emphasized  modeling  shear-producing  dynamic 
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processes  in  the  seasonal  thermocline  in  order  to  aid  in  data 
interpretation  and  to  produce  models  suitable  for  estimating 
shear  levels  under  a  variety  of  oceanographic  and  meteoro¬ 
logical  conditions  (Rubenstein,  1981,  1982).  A  data  analysis 
effort,  described  in  this  report,  presents  statistical  de¬ 
scriptions  of  shear  from  an  array  of  three  current  meter 
moorings . 


Very  few  data  sets  exist  which  address  the  time- 
variability  of  upper  ocean  vertical  shear  on  small  scales. 
Two  recent  sets  are  the  MILE  moorings  and  the  JASIN  moorings. 
Another  set  is  the  one  used  here  which  we  have  chosen  because 
it  includes  three  moorings  separated  by  several  tens  of  kilo¬ 
meters,  and  contains  relatively  long  time  series.  Although 
the  smallest  vertical  separation  is  20  meters,  the  data  are 
particularly  appropriate  for  characterizing  the  near-inertial 
components  of  shear,  and  provide  horizontal  information  that 
the  other  data  sets  do  not. 

This  technical  report  describes  analyses  of  verti¬ 
cal  shear  time-series  using  data  taken  by  the  Woods  Hole 

Oceanographic  Institution  between  27  June  and  18  August  1970 
in  the  vicinity  of  Site  D  (39°10'N,  70°W) .  Three  similar 

surface  moorings  (numbered  338,  339,  and  340),  spaced 

50-70  km  apart,  were  deployed  in  the  geographical  locations 
shown  in  Figure  1.1.  On  each  mooring  were  a  wind  recorder 

and  current  meters  at  12,  32,  52,  and  72  m  depths.  Figure 

1.2  shows  the  distribution  of  data  in  time  for  each  current 
meter.  The  time-series  from  current  meter  3383  was  somewhat 
short,  and  3384  returned  no  data.  All  the  other  time  series 
were  nearly  perfect.  A  schematic  diagram  of  the  mooring 

configuration  is  shown  in  Figure  1.3,  and  Figure  1.4  shows 
the  corresponding  stratification  profiles  presented  by 
Pollard . 
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LOC&  T/ON 


MOORING 

236  39*34. fN,  69°55.5'W 

339  39*07.6  'N,  70*02.3' W 

340  29*07.5'N,  70*35.2’ 


Figure  1.2:  Distribution  of  data  during  Summer,  1970. 

The  first  three  digits  of  the  numbers  on  the 
right  are  the  mooring  numbers.  The  fourth  digit 
indicates  the  sequential  position  of  the 
instrument  on  the  mooring  counting  from  the  top 
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Figure  1.4:  Brunt-Vaisal a  frequency  at  the  three  moorings 
on  9  and  23/24  July  and  16^17  August  as  func¬ 
tions  of  depth  in  the  top  100  m.  The  profiles 
were  computed  using  averaged  values  of  density 
at  approximately  5  m  intervals,  so  details  of 
the  structure  particularly  in  the  top  15  m, 
have  been  smoothed.  From  Pollard  (1980). 
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Pollard  (.1980)  examined  the  velocity  records  Irorr. 
these  data  in  terms  of  the  near-inertial  frequency  oscil¬ 
lations.  He  computed  power  spectra  and  coherences  of 
velocities  m  the  near-inert lal  frequency  range.  He  found 
that  at  12  m  and  32  m,  the  velocities  are  horizontally 
coherent  at  inertial  frequencies  with  suggested  wavelengths 
of  700-1700  km.  At  52  m  and  72  m,  the  records  were  not 
horizontally  coherent,  perhaps  due  to  differences  in 
stratification  among  the  mooring  locations.  Significant 
differences  among  the  12  m  records  were  related  to 
differences  in  local  wind  over  the  50  km  mooring  spacmgs. 
Pollard  found  these  differences  could  be  modeled  using  the 
Pollard  and  Millard  (1970)  model  with  the  local  wind  at  each 
mooring  as  a  forcing  term.  During  periods  of  strong  inertial 
wave  generation,  up  to  two-thirds  of  the  energy  in  the  mixed 
layer  was  at  near  inertial  frequencies.  Finally,  the 
vertical  component  of  the  group  velocity  in  the  top  100  m  was 
downward,  confirming  that  the  flow  of  energy  was  downward 
from  the  surface. 

In  this  report  we  present  calculations  of  power 
spectra,  coherence  spectra,  auto-c orrel at  ions ,  and  cross- 
correlations  from  vertical  shear  time-senes.  Vertical  shear 
is  the  vertical  gradient  of  horizontal  velocity.  For 
purposes  of  our  calculations,  we  approximate  the  shear  by 
computing  velocity  differences  between  pairs  of  vertically 
separated  time-senes  records. 

1.3  THE  GEODYNE  850  CURRENT  METER 

The  Geodyne  850  Current  Meter  measures  current 
speed  using  a  Savonius  rotor,  and  current  direction  using  a 
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freely  moving  directional  vane  and  a  magnetic  compass.  Tbe 
measurements  were  recorded  in  the  "interval"  mode.  lr.  this 
mode,  bursts  of  rapid  samples  were  recorded  every  15  minutes. 
Each  recording  interval  consisted  of  a  burst  of  20  samples  at 
5-second  intervals,  followed  by  a  longer  interval  (about  800 
seconds)  during  which  no  samples  were  taken.  The  currents 
were  recorded  as  speed  and  direction.  As  part  of  the  data 
processing,  they  were  transferred  into  Cartesian  coordinates, 
and  then  vector  averaged  to  obtain  a  single  measurement  every 
1 5  minutes. 

The  response  of  the  Savonius  rotor  is  not  well 
understood.  There  is  evidence  that  the  rotor  accelerates 
faster  than  it  decelerates  (Fofonof  and  Ercan,  1967),  but 
there  is  also  evidence  of  the  opposite  being  true  (see 
response,  Halpern  et  al. ,  1981).  In  either  case,  tbe  time 
constants  of  the  rotor  and  directional  vane  are  not  perfectly 
matched . 

A  typical  mooring  (338)  is  shown  schematically 
in  Figure  1.3.  The  other  two  moorings  were  similar, 
especially  the  upper  four  current  meters  with  which  we  are 
concerned.  The  toroidal  surface  buoys  support  a  metal 
triangular  pyramid  superstructure.  Within  this  super¬ 
structure  is  placed  the  wind  recorder.  In  order  to  prevent 
an  upwind  superstructure  strut  from  creating  spurious 
effects,  a  large  plywood  fin  is  attached  to  one  of  the 
struts.  This  helps  to  orient  the  buoy  so  that  the  wind 
impinges  on  the  wind  sensor  midway  between  two  struts. 

Halpern  et  al.  (1981)  performed  tests  to  determine 
the  effects  of  mooring  motions  on  current  meter  measure¬ 
ments.  He  compared  the  spectra  from  upper  ocean  (above 
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100  ns )  current  meters  supported  by  surface  buoys  with  current 
meters  supported  by  spar  and  subsurface  buoys.  He  defined  a 
quantity  called  the  "maximum  usable  frequency",  the  frequency 
threshold  above  which  the  power  spectra  from  two  current 
meters  (at  the  same  depth  but  attached  to  different  moorings) 
diverge.  In  the  upper  30  m,  the  maximum  usable  frequency  was 
at  least  2  cph  m  all  cases.  In  the  depth  range  79-94  m, 
where  the  currents  were  much  weaker,  the  maximum  usable 
frequency  fell  m  the  range  0.2-2  cph.  As  we  are  interested 
in  shear  near  the  inertial  frequency,  0.053  cph,  this 
limitation  is  of  little  consequence. 


1-9 


r 


Section  2 

STATISTICS  PROCEDURES 


2.1  POWER  SPECTRAL  DENSITIES 

Calculations  of  power  spectral  densities  follow  the 
procedure  outlined  below. 

1.  Select  a  length  of  record  with  a  number  of  data  points  M 
equal  to  a  power  of  two. 

2.  The  data  record  is  multiplied  by  a  Hanning  window 
function 

WT/2(t)  »  |  A(1  +  cos  2nt/T), 

where  t  is  time,  T  is  the  total  length  of  the  record,  and  the 

scale  factor  A  -  /873  in  sures  that  the  variance  remains 
unchanged.  The  corresponding  filter  W(f)  is  given  by  the 
Fourier  transform  of  wT^2(t); 

f(f)  »  S1P**T  (l-f2T2)-1. 

Figure  2.1  shows  the  shape  of  W( f ) ,  for  the  particular  case 
T  -  1. 

3.  The  fast  Fourier  transform  of  the  u-component  is  computed 
as  follows: 

Dk  .  FFT(um)  .  “f  u  exp  {-i  ^$E) 
m=o 
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Figure  2.1  Frequency  response  W(f) 

for  a  Hanning  window  function. 
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*  kth  spectral  coefficient 
Up,  *  mtb  u-component  of  velocity 

and  the  v-component  is  computed  similarly. 

4.  Raw  spectral  estimates  are  computed; 


corresponding  to  frequencies 

fk  *  AtM  ;  k  =  °’ 1 ’ * ‘ * ’“Z2-1 • 

5.  Band  averaging  is  performed,  providing  50  smoothed 
spectral  estimates.  The  number  of  raw  estimates  used  to 
compute  the  smoothed  spectrum  varies  from  2  in  the  lowest 
bins  up  to  347  m  the  highest  bin.  The  distribution  of 

-smoothed  estimates  is  designed  to  keep  the  values  of  log  fk 
nearly  evenly  spaced,  and  to  resolve  the  inertial  and  semi¬ 
diurnal  frequencies.  In  the  bin  containing  the  inertial 
frequencies,  7  raw  estimates  are  averaged. 

6.  The  smoothed  spectral  estimates  are  mutilplied  by 
frequency  f,  and  plotted  on  a  linear  scale  versus  log  f. 
This  format,  often  called  a  variance-preserving  riot, 
provides  constant  variance  per  unit  area  under  the  curve. 
Frequency  bands  which  contribute  significantly  to  the  total 
variance  are  easily  recognized  by  this  method. 
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2.2  VECTOR  COHERENCE  SPECTRAL  FUNCTIONS 

®  We  follow  the  procedure  outlined  by  Gonella  (1972) 

for  computing  the  vector  coherence  functions.  This  procedure 
is  equivalent  to  that  described  below. 

#  1.  Begin  with  two  discrete  vector  time  series,  in  complex 

form ; 


-  W  +  ivl(t 

W  ’  u2(tm>  +  iv2(t 


m 

m 


), 

); 


ro  *  0,1,..., L- 1 . 


The  integer  L  is  chosen  to  be  an  integer  (I)  multiple  of  J, 
where  J  is  a  power  of  two; 


L  *  IJ. 


Divide  each  of  the  two  time  series  into  21-1  short  series, 
overlapping  one  another  by  50%,  and  each  of  length  J. 

2.  Multiply  each  of  the  short  series  by  a  Hanning  window 
function,  decribed  in  the  previous  section. 

3.  The  two-sided  Fourier  transform  (the  transform  is  not 
symmetric  with  respect  to  zero  frequency)  is  computed  for 
each  short  complex  series; 

Ii<v  -  "V  «p(-i  a^B)  , 

m*o 

k  -  -M/2,  .  .  .  ,  ( M/2 )-l , 
and  similarly  for  Ug( fk) • 
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4.  The-  auto,  coincident,  and  quadrature  spectral  density 
functions  P  and  Q  are  computed: 


P11(V  "  <W 


P22(fk)  -  <W  ^2(1k)> 
P12(fk)  -  xQ12Uk)  =  <Ul(fk)  U*(ik)> 


where  the  asterisk  denotes  a  complex  conjugate,  and  the 
brackets  denote  ensemble  averages  over  the  21  -  1  ensembles. 
This  procedure  provides  estimates  with  21  degrees  of  freedom. 


5.  Compute  the  coherence  amplitude  and  phase,  defined  by 


C<V 

<Kfk) 


12 


+  Q 


12 


P  P 

li  22 


1/2 


* 


arctan(Q12/P12) . 


for  frequencies  fk  •  k/(NAt)»  for  k  = -M/ 2 , . . . , M/ 2- 1 . 
Negative  frequencies  indicate  clockwise  rotation  and  positive 
frequencies  indicate  counterclockwise  rotation. 


2.3  NORMALIZED  CORRELATION  FUNCTIONS 


Normalized  auto-correlat ion  and  cross-correlation 
functions  are  computed  following  the  procedure  outlined 
below. 


1.  A  maximum  lag  g  is  chosen. 


t 


* 
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2.  Normalized  correlation  is  computed: 


C  (  r  Li )  = 

xy 


M-g 

m=l 


l!W  -  * 


“-je  _  2  “-e 
lxm  -  *>  r 

m=l  m=l 


r  =  -g 


>g< 


The  overbars  denote  mean  values. 


» 
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Section  3 
POWER  SPECTRA 


In  this  section  we  discuss  power  spectra  of  the 
velocity  and  shear  (velocity  difference)  time  series.  The 
method  of  computing  power  spectra  is  described  in  Section  2. 
Our  purpose  is  to  examine  the  relative  contributions  by 
individual  frequency  bands.  It  is  convenient  to  categorize 
these  bands  according  to  the  physical  processes  with  which 
they  are  associated.  We  divide  the  frequency  domain  into 
four  bands:  low  frequency,  near-inertial,  semidiurnal  and 
internal  wave.  In  Section  3.1  we  present  a  brief  background 
review  of  the  physical  processes  that  generate  kinetic  energy 
in  these  frequency  bands.  In  Section  3.2  we  discuss  the 
velocity  spectra  computed  by  Pollard  (1980).  In  Section  3.3 
we  present  and  discuss  spectra  of  the  shear  between  various 
pairs  of  current  meters. 

3.1  FORCING  MECHANISMS 

3.1.1  Near-Inertial  Oscillations 


We  shall  see  that  the  near-inertial  frequency 
band  is  generally  a  strong  contributor  to  the  total  shear 
variance  in  upper  ocean  power  spectra.  The  inertial 
frequency  f  is  defined  as  twice  the  magnitude  of  the  pro¬ 
jection  of  the  earth's  angular  rotation  vector  on  the  local 
vertical  unit  vector  (i.e.,  f  -  2QsinX,  where  0  is  the 
earth's  angular  rotation  rate,  and  x  is  the  latitude).  At 
the  latitude  of  Pollard's  current  measurements  (X  *  40°), 
the  inertial  frequency  is  f  *  9.3  x  10~5rad/sec  or  .053 

cycles/hr . 
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The  inertial  frequency  is  the  natural  frequency  of 
oscillation  in  a  rotating  unstratified  fluid.  When  a  sudden 
change  in  the  wind  stress  is  applied  over  the  ocean,  the 
horizontal  currents  in  the  upper  ocean  are  deflected  in  the 
clockwise  sense  (in  the  Northern  Hemisphere)  by  the  Coriolis 
force.  The  velocity  vector  begins  rotating  in  the  clockwise 
direction.  These  oscillations  are  often  observed  in  the 
upper  ocean,  for  several  days  after  the  passage  of  an 
atmospheric  front. 

3.1.2  Semi-Diurnal  Oscillations 


The  oscillations  observed  at  the  semi-diurnal 
frequency  are  tidal  currents.  Regal  and  Wunsch  (1973) 
describe  observations  of  internal  tides  at  several  upper 
ocean  stations,  including  Site  D.  At  this  site,  most  of  the 
semidiurnal  energy  is  in  baroclinic  modes,  and  is  classified 
as  the  lunar  (M2)  constituent,  with  a  period  of  12.42 
hours.  Regal  and  Wunsch  explain  that  the  baroclinic  tidal 
component  is  predominant  here  because  of  the  interaction  of 
tidal  motions  with  the  geometry  of  the  nearby  edge  of  the 
continental  shelf.  Internal  tides  are  generated  at  the  sharp 
break  in  the  bottom  slope,  just  to  the  north  of  Site  D.  As 
shown  in  Figure  3.1,  internal  wave  energy  propagates  along 
characteristic  lines  given  by  z+cy  *  constant,  where  z 
is  vertical,  y  is  horizontal,  and  c2  =  (u>2  -  f 2 )  / 
(N2(z)  -  u2).  One  of  the  internal  wave  characteristics, 

represented  by  the  solid  line  labeled  A,  lies  nearly  tangent 
to  the  continental  slope,  and  reaches  Site  D  after  a  single 
bounce.  Regal  and  Wunsch  suggest,  without  detailed  explana¬ 
tion,  that  internal  wave  energy  is  focused  at  the  surface 
near  the  site.  Other  characteristics  also  arrive  near  Site 
D,  but  the  focusing  effect  is  not  as  pronounced. 
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Figure  3.1:  Profile  of  relief  along  70°W  (see  Figure  1.1 
for  local  isobaths)  together  with  selected 
semidiurnal  period  internal  wave  character¬ 
istics  passing  near  site  D.  Note  the  solid 
characteristic  line  (labeled  A)  which  is 
nearly  tangent  to  the  continental  slope  and 
reaches  site  D  in  a  single  bounce.  (Regal 
and  Vunsch .  1973) 
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3.1.3 


Low  Freouencv  Motions 


We  have  partitioned  the  low-  frequency  band  as  uK  f 
because  it  represents  motions  caused  by  several  different 
physical  processes  excluding  internal  waves.  Internal  waves 
cannot  propagate  freely  at  frequencies  less  than  the  local 
inertial  frequency.  Low  frequency  motions  have  a  variety  of 
generating  mechanisms.  Among  these  mechanisms  are  slow 
changes  in  the  large-scale  wind  stress  and  changes  in  the 
ocean's  baroclimc  density  structure.  We  will  not  further 
describe  these  processes  in  this  report.  Our  shear  spectra 
generally  do  not  show  much  energy  in  the  low  frequency  bands. 

3.1.4  Internal  Wave  Band 


A  large  number  of  mechanisms  have  been  advanced 
for  the  generation  of  internal  waves.  We  touch  upon  just  a 
few  of  the  most  likely  mechanisms  here.  The  variability  of 
internal  -  wave  energy  in  the  upper  ocean  suggests  that  some 
mode  of  surface  forcing  may  be  important.  Surface  forcing  by 
traveling  fluctuations  of  divergent  wind  stress  induces 
vertical  motions  at  the  bottom  of  the  surface  boundary  layer. 
Traveling  fluctuations  in  buoyancy  flux  give  rise  to  pressure 
fluctuations  which  also  induce  vertical  motions  at  the  base 
of  the  surface  layer.  These  motions  generate  internal  waves 
which  are  free  to  propagate  into  the  interior  (Thorpe,  1975). 
Internal  waves  may  be  induced  by  currents  flowing  over  bottom 
topography.  We  have  already  discussed,  in  Section  3.1.2,  the 
possibility  of  internal  waves  being  generated  by  tidal 
motions  through  the  latter  mechanism. 

In  this  report  we  propose  yet  another  possible 
mechanism  for  the  generation  of  some  portion  of  the  kinetic 


energy  m  me  internal  wave  band.  Under  certain  conditions, 
strong  inertial  oscillations  may  interact  nonlmearly  with 
the  strong  tidal  component.  Energy  at  the  difference 
frequency,  being  smaller  than  f,  cannot  propagate,  but 
energy  at  the  sum  of  the  inertial  and  semidiurnal  fre¬ 
quencies  can  propagate  as  a  free  internal  wave  into  the 
stratified  interior.  Our  spectra  show  evidence  of  a 
component  at  this  frequency.  Krauss  (1981)  observed  a 
similar  spectral  component  in  the  Baltic  Sea  at  twice  the 

inertial  frequency. 

3.2  VELOCITY  SPECTRA 

Figure  3.2  shows  the  velocity  spectra  computed  by 
Pollard  (1980).  Two  peaks  at  the  inertial  and  semidiurnal 

frequencies  are  the  most  prominent  features  in  these 

spectra.  At  the  two  shallowest  depths,  12  and  32  m,  the 
inertial  peak  is  stronger  relative  to  the  semidiurnal  peak. 
At  '52  and  72  m,  the  opposite  is  the  case.  Note  that  the 

vertical  scales  of  the  plots  are  different  at  each  depth. 

The  semidiurnal  peak  is  most  energetic  at  all 
depths  at  moorings  339  and  340.  As  discussed  m  Section 

3.1.2,  this  may  be  due  to  their  position  relative  to  the 

continental  slope.  Conditions  are  such  that  moorings  339  and 
340  are  in  the  most  favorable  location  to  receive  the 
benefits  of  the  internal  tides  focusing  effect. 

In  the  mixed  layer,  at  12  m,  the  inertial  peaks 

recorded  at  the  three  moorings  were  of  comparable  magnitude. 
At  32,  52,  and  72  m,  mooring  340  shows  the  strongest  in¬ 
ertial  oscillations.  Pollard  (1980)  showed  that  two  fronts 
passed  each  of  the  current  meter  arrays  on  11  August,  and 

that  during  this  period  current  meter  3383  was  not  in 
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Figure  3.2:  Spectra  of  all  records  in  the  top  100  m  using  1024  h  of 
data  between  28  June  and  10  August.  The  abscissa  is 
log-frequency  on  the  lower  scale,  log-period  on  the  upper 
»  scale  of  each  plot.  The  ordinate  is  the  energy  density, 

scaled  by  the  frequency,  plotted  linearly.  Note  that  the 
energy  scale  is  different  at  different  depths,  but  the  same 
for  all  records  at  a  fixed  depth.  The  area  under  this  plot 
in  any  given  frequency  band  is  proportional  to  the  energy 
in  that  band.  Dominant  peaks  are  inertial  (18.78  h)  and 
•  tidal  (12.42  h).  From  Pollard  (1980). 
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operation.  The  time  difference  between  arrival  of  the  fronts 
at  340,  338,  anc  339  were  7.  12,  and  16  hours,  respec t ivel y . 
As  a  result,  the  inertial  amplitudes  after  the  passage  of  the 
fronts  were  observed  to  be  much  stronger  at  340  than  at  339. 
Ax  339,  energy  input  by  the  wind  during  the  first  half  of  an 
inertial  period  is  subtracted  during  the  second  half.  Thus, 
the  energy  at  the  inertial  frequency  appears  to  be  quite 
sensitive  to  small  differences  in  the  wind  stress. 

The  relative  contributions  by  the  low  frequency 
and  internal  wave  bands  are  practically  negligible  at  12  and 
32  m.  At  52  and  72  m,  the  relative  contributions  in  these 
frequency  bands  are  somewhat  greater. 

3.3  SHEAR  SPECTRA 

3.3.1  Resul ts 

Power  spectral  density  was  computed  for  the 
velocity  differences  between  pairs  of  vertically  adjacent 
current  meters.  Since  the  interval  between  data  samples 
was  15  minutes,  the  Nyquist  frequency  was  2  cycles  per 
hour.  Since  the  fast  Fourier  transform  algorithm  requires 
■the  record  length  to  be  a  power  of  two,  the  record  length  was 
chosen  to  contain  4096  samples.  The  data  were  not  filtered 
prior  to  the  computation  of  spectral  estimates.  Table  3.1 
lists  the  current  meter  pairs  for  which  spectra  were 
computed,  along  with  their  respective  depths.  The  spectral 
density  functions  are  displayed  m  variance-preserving  format 
in  Figures  3.3  through  3.9.  Velocity  differences  between 
non-adjacent  current  meters  were  not  analyzed  because  these 
differences  are  dominated  by  the  relatively  large  velocities 
at  the  upper  current  meter.  Analysis  of  these  velocity 
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TABLE  3.1 


POWER  SPECTRA  OF  VELOCITY  DIFFERENCES 
(VERTICAL  SHEAR) 


Fig.  Current  Meters  Depths  (m) 


3.3 

3393-3394 

12-32 

3.4 

3402-3403 

12-32 

3 . 5 

3394-3395 

32-52 

3.6 

3403-3404 

32-52 

3.7 

3385-3386 

52-72 

3.8 

3395-3396 

52-72 

3.9 

3404-3406 

52-72 
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differences  would  Therefore  be  nearly  equivalent  to  the 
analysis  of  individual  current  meter  records  previously  done 
by  Pollard  (1980). 

The  spectra  in  Figures  3.3  -  3.9  are  for  horizon¬ 
tal  velocity  differences  over  Az  =  20  m  vertical  separation. 
We  can  convert  these  to  shear  spectra  simply  by  dividing  by 
hz2.  Grabowski  and  Newman  (1981)  used  a  vertical  wave- 
number  shear  spectral  model  based  on  that  compiled  by 
Gargett  et  al.  (1981),  to  determine  the  dependence  of  mean 
square  shear  on  Lz.  Grabowski  and  Newman  found  that  mean 
square  shear  varies  as  Az-1.  Therefore,  much  of  the 
small  vertical-scale  shear  variance  is  missed  due  to  the 
large  (20  m)  vertical  separation  of  the  current  meters. 
Extending  the  results  m  Table  1  of  Grabowski  and  Newman,  the 
shear  variance  observed  with  a  vertical  sample  interval  Lz  = 
20  m  is  54%  of  the  variance  observed  with  Lz  *  10  m  and  14% 
of  the  variance  observed  with  Az  =  2  m. 

3.3.2  Inertial  and  Tidal  Shear 

Sharp  peaks  are  evident  m  all  of  the  spectra  at 
the  inertial  frequency,  .053  cpb.  The  peaks  are  strongest 
at  the  shallowest  depth  pair  (12  -  32  m) ,  and  decrease  in 
magnitude  with  depth.  As  was  the  case  for  the  velocity 
spectra,  the  inertial  band  energy  varies  from  one  mooring  to 
the  next,  because  the  shear  response  is  sensitive  to  small 
differences  in  the  wind  stress. 

The  inertial  peak  in  each  (20  m  velocity  differ¬ 
ence)  shear  spectrum  is  comparable  in  magnitude  to  the  peak 
in  the  velocity  spectrum  of  the  corresponding  upper  current 
meter.  The  reason  is  that  the  inertial  oscillations  decay 
rapidly  with  depth.  The  20  m  velocity  differences  are 
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therefore  dominated  by  the  velocities  recorded  by  the  upper 
current  meters. 


Distinct  peaks  at  the  semidiurnal  frequency 
are  present  in  most  of  the  shear  spectra,  but  only  in  Figure 
3.4  is  the  peak  comparable  to  the  inertial  peak.  This 
result  is  in  contrast  with  the  velocity  spectra,  which  show 
strong  peaks  of  magnitude  comparable  to  the  inertial  peaks 
for  about  half  of  the  records.  The  implication  is  that, 
even  if  the  tidal  oscillations  are  baroclinic,  their  verti¬ 
cal  wavelengths  are  generally  much  greater  than  20  m. 
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Internal  Wave  Band 


The  spectral  levels  of  vertical  shear  increase 
with  frequency  from  0.5  -  2.0  cph  iD  the  internal  wave  band. 
In  Pollard's  velocity  spectra  (Figure  3.2)  this  increase  is 
not  apparent  because  the  spectra  were  terminated  at  0.5  cph. 
This ‘  increase  may  .  be  a.  contaminating  effect  of  mooring 
motion,  an  aliasing  effect,  or  some  combination  of  both. 

Halpern  et  al .  (1S81)  discuss  intercomparison  tests  between 

current  measurements  in  the  upper  ocean  made  with  different 
types  of  current  meters  and  moorings.  They  find  that  VACM 
(Vector  Averaging  Current  Meter)  measurements  are  more 
susceptible  to  spurious  mooring  motions  in  the  seasonal 
thermocline  (below  50  m)  than  in  the  mixed  layer  where  the 
current  magnitudes  are  stronger.  In  comparisons  between 
VACM  measurements  taken  from  moorings  attached  to  surface 
and  subsurface  buoys,  the  surface  buoy  data  are  more  noisy. 
Halpern  et  al  .  define  a  "maximum  usable  frequency"  as  the 

frequency  above  which  two  compared  spectra  begin  to  diverge. 
In  the  mixed  layer,  the  maximum  usable  frequency  is  about  4 
cycles  per  hour.  In  the  seasonal  thermocline  (79  m)  the 

maximum  usable  frequency  ranges  from  about  0.2  to  2  cycles 
per  hour. 


In  the  present  data  set,  the  current  teters 
recorded  data  at  the  rate  of  4  samples  per  hour.  It  is 
possible  that  the  data  were  under-samplec:  energy  at  higher 
frequencies  may  be  aliased,  or  folded  into  the  computed 
spectra.  We  cannot  judge  the  degree  of  aliasing  from  these 
spectra  alone.  However,  we  can  show  an  example  of  how- 
aliasing  affects  a  model  spectrum.  Figure  3.10  shows  a 
portion  of  a  spectrum  of  horizontal  velocity,  based  on  the 
Garrett  and  Munk  model  (Munk,  1981).  The  spectrum  is  shown 
in  variance  preserving  f orm--variance  in  a  given  frequency 
band  is  proportional  to  the  area  under  the  curve.  It  is 
based  upon  an  inertial  frequency  f  =  0.053  cph  and  a  local 
Vaisala  frequency  N  =  4  cph.  If  the  current  meter  sampling 
rate  is  4  samples  per  hour,  then  the  Nyquist  frequency  is 
2  cph.  Power  in  the  frequency  band  between  2  and  4  cph  is 
aliased  into  the  computed  spectrum.  The  aliased  spectrum 
overestimates  the  spectral  magnitude,  especially  between  1 
and  2  cph. 


3.3.4  Dissipation  of  Inertial  Oscillations 

One  of  the  biggest  problems  in  our  understanding 
of  inertial  oscillations  is  the  mechanism  by  which  they  are 
dissipated  in  the  upper  ocean.  Several  mechanisms  have  been 
suggested,  but  none  of  these  prove  to  be  wholly  satisfactory. 

Pollard  (1980)  estimates  that  in  the  mixed  layer, 
the  inertial  oscillations  during  the  period  30  June  -  5  July 
decreased  in  amplitude  from  40  cm  s-1  to  20  cm  s_1.  This 
corresponds  to  an  energy  loss  rate  of  4  ergs  cm-2  s-l.  He 
derives  a  simple  estimate  of  the  energy  loss  rate  due  to 
simple  dispersion  of  near-inertial  internal  waves,  approxi¬ 
mately  0.5-1. 5  erg  cm-2  s-1.  This  discrepancy  shows  that 
dispersion  alone  cannot  account  for  the  energy  loss. 
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.05  .10  .  25  .  50  1.00  2.00  5.00 

Cycles  Per  Hour 

Figure  3.10:  Garrett  and  Munk  horizontal  velocity  spectral  model GM79) 
frequency  dependence  Fu(u>)  «  (  w2  +  f2)ur3(u;2  -  f2)“-'- 
in  variance  preserving  form.  Inertial  frequency  is  1  = 
.053  cph,  Vaisala  frequency  is  N  =  4  cph.  Also,  if  data 
are  under-sampled  at  a  rate  of  4  samples  per  hour,  the 
spectrum  is  aliased  as  shown. 
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The  generation  time  scale  of  inertial  oscilla¬ 
tions  (~1  cay)  is  much  shorter  than  the  dissipation  time 
scale  (~5  days).  This  fact  helps  to  explain  why  the  oscill¬ 
ations  are  concentrated  in  the  surface  layer.  Widely  spacec 
wind  events  quickly  stir  up  inertial  motions,  which  take 
many  days  to  dissipate. 

Other  mechanisms  for  the  dissipation  of  inertial 
energy  have  been  advanced.  Pollard  et  al .  (1973)  suggest 

the  possibility  that  the  Kelvin -Helmholtz  instability  at  the 
mixed  layer  base  ma-  remove  energy  from  the  inertial  oscill¬ 
ations.  This  suggestion  is  supported  by  observations  and 
numerical  computations  of  small  Richardson  number  at  the 
mixed  layer  base  by  Krauss  (1981).  Bell  (1978)  suggests 
that  turbulence-induced  undulations  in  the  mixed  layer- 
thermocline  interface  may  be  advected  by  inertial 
oscillations.  These  undulations  are  forced  to  do  work 
against  the  underlying  stratified  fluid,  and  may  generate 
radiation-damping  internal  waves. 

Our  shear  spectra  suggest  that  another  mechanism 
may  also  be  at  work.  There  are  small  elevations  in  the 
spectral  levels  at  the  frequency  ~0.13  cpb ,  the  sum  of  the 
inertial  and  semidiurnal  frequencies.  Their  consistent 
appearance  in  most  of  the  spectra  is  suggestive  of  a  non¬ 
linear  resonant  interaction  taking  place.  If  a  nonlinear 
interaction  is  sufficiently  strong,  then  energy  at  the  sum 
and  difference  frequencies  f  _+  u>g ,  (where  ^  is  t'3e  sem*“ 
diurnal  frequency)  is  generated.  Energy  at  the  difference 
frequency,  being  smaller  than  f,  is  not  free  to  propagate. 
Energy  at  the  sum  frequency  is  free  to  propagate  into  the 
stratified  interior.  As  the  angle  of  propagation  of  this 
internal  wave  component  is  steeper  than  the  near-inertial 
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component,  its  vertical  group  velocity  is  greater  anc  hence, 
its  energy  dispersion  rate  may  also  be  greater.  Observa¬ 
tions  by  Krauss  (1981)  in  the  Baltic  Sea  show  spectral  peaks 
at  2 f.  A  similar  resonant  nonlinear  interaction  may  be 
occurring  there,  and  may  be  further  evidence  for  such  a 
radiation  damping  mechanism. 


Section  4 
COHERENCE  SPECTRA 

Coherence  spectral  functions  were  computed  between 
pairs  of  vertical  shear  measurements  (velocity  differences), 
separated  horizontally  and  vertically.  Coherence  and  phase 
spectra  were  also  computed  between  wind  stress  and  vertical 
shear  time  series.  Section  2  described  the  procedure  used  to 
compute  vector  coherence  and  phase. 

4.1  HORIZONTAL  COHERENCES 

Following  the  notation  of  Section  2,  all  coherence 
computations  were  based  on  records  of  length,  L  *  4608.  In 
the  horizontal  and  slant  coherences,  and  in  the  wind  stress/ 
shear  coherence  calculations,  most  of  the  significant 
coherence  was  in  the  frequency  range  (-.25,. 25)  cph .  Here 
negative  frequencies  indicate  clockwise  rotation,  positive 
frequencies  indicate  counterclockwise  rotation.  In  order  to 
provide  sufficient  resolution  in  this  frequency  range,  M  was 
chosen  to  be  9,  giving  17  overlapping  series,  each  of  length 
J  *  512.  This  choice  allowed  2M  =  18  degrees  of  freedom. 
The  90%  confidence  level  was  0.53. 

Horizontal  coherence  amplitude  and  phase  functions 
were  computed  between  vertical  shear  (velocity  difference) 
time  series.  Table  4.1  conveniently  lists  the  figures,  with 
their  respective  current  meter  numbers,  depths,  and  approxi¬ 
mate  horizontal  separations.  The  coherence  spectra  are 
displayed  in  Figures  4.1  through  4.5.  The  frequency  range  of 
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TA3LE  4.1 

SPECTRA  BETWEEN  HORIZONTALLY 
VERTICAL  SHEAR  MEASUREMENTS 


Approximate 

Horizontal 


• 

Fig. 

Current  Meters 

Depth  (m) 

Separat; 

4.1 

3393-3394/3402-3403 

12-32 

47 

4.2 

3394-3395/3403-3404 

32-52 

47 

• 

4.3 

3395-3396/3404-3406 

52-72 

47 

4.3 

3395-3396/3385-3386 

52-72 

51 

4.5 

3404-3406/3385-3386 

52-72 

69 
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Figure  4.5  Horizontally  Separated  Vertical  Shear  Coherence 
Amplitude  (Top)  and  Phase  (Bottom), 

Between  3404-3406  and  3385-3386. 
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display  is  between  -0.25  to  0.25  cph.  In  Figures  4.1  and 
4.2,  the  coherence  amplitude  is  strong  in  the  frequency  range 
(-f,  0).  This  result  is  similar  to  Gonella's  (1972)  finding 
that  the  coherence  between  horizontally  separated  near- 
surface  current  measurements  was  strong  in  this  range.  For 
deeper  shear  measurements,  the  coherence  amplitude  (Figures 
4.3  through  4.5)  is  also  significant  at  some  frequencies  in 
this  range,  but  it  is  not  as  strong  as  for  the  near-surface 
measurements . 

4.2  VERTICAL,  COHERENCES 

Vertical  rotary  coherence  spectra  of  shear  time 
series  were  computed,  but  showed  rather  anomalous  results. 
For  example,  some  rotary  spectra  of  vertically  adjacent  time 
series  exhibit  a  minimum  in  coherence  at  the  inertial  fre¬ 
quency.  These  calculations  and  their  interpretation  require 
further  study. 

4.3  WIND  STRESS/ SHEAR  COHERENCES 

Coherence  amplitude  and  phase  functions  were 
computed  between  wind  stress  and  vertical  shear  time  series. 
Table  4.2  lists  the  wind  recorders  and  current  meters  used  in 
the  coherence  calculations. 

Wind  stress  _x  was  computed  from  the  drag  law 

J  *  PCU  j  U  j , 

where  p  is  the  air  density,  and  U  is  the  mean  wind  velocity 
vector,  with  magnitude  jUj  .  The  stress  coefficient  C  was 


TABLE  4.2 

COHERENCE  SPECTRA  BETWEEN  WIND  STRESS 
AND  VERTICAL  SHEAR  MEASUREMENTS 


Fifi. 

Wind  Recorder/Current  Meters 

Depth  (m) 

4.6 

3391/3393-3394 

12-32 

4.7 

3401/3402-3403 

12-32 

4.8 

3391/3394-3395 

32-52 

4.9 

3401/3403-3404 

32-52 

4.10 

3381/3385-3386 

52-72 

4.11 

3391/3395-3396 

52-72 

4.12 

3401/3404-3406 

52-72 
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computed  using  the  relation  derived  empirically  by  Amorocho 
and  DeVries  (1980): 


C  =  0.0015  [l  +  exp(-W  1~53'2,5  )  l'1  +  0.00104 
4.4  COHERENCE  ANALYSIS 

We  center  our  discussion  of  coherence  on  the 
coherence  spectra  of  wind  stress  and  vertical  shear  displayed 
in  Figures  4.6  -  4.12.  The  shallowest  shear  records  (12  - 

32m)  are  coherent  with  the  wind  stress  in  the  frequency  band 
0  >  u>  >  -f  (see  Figures  4.6  and  4.7).  In  the  vicinity  of 
the  inertial  frequency,  the  phase  shifts  by  approximately 
180°.  This  phase  shift  if  expected  because  the  surface 
layer  is  resonant  at  the  inertial  frequency.  In  the  next 
pair  of  coherence  spectra  (Figures  4.8  and  4.9)  between 
wind  stress  and  shear  at  32  -  52m,  the  coherence  amplitude  is 
not  significant  at  the  inertial  (f  =  -0.053  cph).  Figure 
4.  8  shows  no  significant  coherence  at  the  lowest  frequencies 
|u)J  <  f.  Figure  4.9  shows  a  slight  degree  of  coherence  in 
the  frequency  band  |w|  <  f/2.  Of  the  three  coherence  spectra 
at  the  deepest  levels  (Figures  4.10  -  4.12;  52  -  72m),  none 
show  coherence  peaks  at  the  inertial  frequency. 

The  coherence  features  discussed  thus  far  are  very 
similar  to  those  found  by  Weller  (1981).  Weller  observed  the 
upper-ocean  current  by  lowering  vector  measuring  current 
meters  from  FLIP.  Whereas  in  our  data  regime  the  mixed  layer 
depth  varied  from  10  to  20m  (see  Figure  1.4),  Weller's 
observations  were  taken  in  a  regime  where  the  mixed  layer 
depth  was  about  50  m.  Figure  4.  13  shows  coherence  between 
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Figure  4.8  (Left):  Between  3391  and  3394-3395. 
Figure  4  .9  (Right):  Between  3401  and  3403-3404. 
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FREG  (CPH)  «!0"  FREQ  ( CPH )  «10" 

Wind  Stress/Vertical  Shear  Coherence  Amplitudes  (Top)  and  Phases  (Bottom) 
Figure  4,lo  ("Left):  Between  3381  and  3385-3386. 

Figure  4.11  (Right):  Between  3391  and  3395-3396. 
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Figure  4.12  Wind  Stress/Vertical  Shear  Coherence 
Amplitudes  (Top)  and  Phases  (Bottom) 
Between  3401  and  3404-3406. 


wind  stress  and  shear,  as  a  function  of  depth  and  frequency. 
In  this  figure,  the  shear  in  the  mixed  layer  is  highly 
coherent  in  the  frequency  band  0  >  >  -f.  The  coherence 

extends  deepest  at  the  lowest  frequencies,  and  at  the 
inertial  frequency.  At  the  inertial  frequency,  there  is  no 
significant  coherence  between  50  and  80  m,  but  below  80  m, 
there  is  again  significant  coherence.  Weller  suggests  down¬ 
ward  propagating  near-inertial  waves,  caused  by  local  wind- 
induced  Ekman  suction  as  a  possible  explanation  for  this 
feature.  It  may  also  be  possible  that  fluid  in  the  region  of 
the  sharp  density  interface  just  below  the  surface  layer  is 
vertically  advected  past  the  current  meters.  This  could 
result  in  shear  records  with  discontinuities  which  tend  to 
reduce  the  coherence  with  wind  stress.  In  Rubenstein 
(1981a, b)  we  further  discuss  this  mechanism,  and  give 
quantitative  descriptions  of  internal  wave  propagation. 

In  addition  to  high  coherence  in  the  range 
-f  <  u  <  0,  two  other  interesting  features  show  up  in  Figures 
4.  6  and  4.7  .  In  Figure  4.  9  there  is  a  coherence  peak 
which  just  reaches  the  significance  level  at  u  «  2f  =  -.105 
cph.  In  Figure  4.7  there  is  a  strong  coherence  peak  at 
u  =  f  +  =  -.138  cph  (where  denotes  the  semidiurnal 

frequency).  We  discussed  in  Section  3  the  evidence  for 
nonlinear  interactions  between  inertial  oscillations  and 
semidiurnal  tides.  These  nonlinear  interactions  showed  up  in 
the  shear  power  spectra  as  small  peaks  at  o>  =  f  +  Wg  (see 
Section  3.3.4).  We  are  probably  seeing  the  coherence  between 
some  combination  of  the  following  nonlinear  effects: 


r 


r*eoucNcr 


FRCQUCNCY 


Figure  4.13  Coherence  (top)  and  phase  (bottom)  between 
the  wind  stress  and  the  shear.  Positive 
angles  indicate  shear  is  to  the  right  of  the 
<•  wind  stress.  From  Weller  (1981). 
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1.  There  are  usually  strong  shears  associated  with 
the  sharp  density  interface  at  the  base  of  the 
mixed  layer.  These  shears  are  not  likely  to  be 
linear  in  depth.  A  pair  of  current  meters 
located  in  the  vicinity  responds  nonlinearly  to 
small  vertical  displacements  of  the  interface. 

2.  We  use  an  imperfect  parameterization  in  cal¬ 
culating  the  wind  stress.  The  deviation  of 
this  parameterization  from  the  true  wind  stress 
is  a  nonlinear  function  of  the  wind  speed. 

3.  The  motion  of  the  surface  buoy  responds  in  a 
nonlinear  way  to  changes  in  the  direction  of 
the  wind  and  surface  current.  Information 
about  the  motion  of  the  buoy  may  be  transmitted 
through  the  mooring  line  to  the  individual 
current  meters.  These  effects  may  be  compli¬ 
cated  even  further  by  motions  caused  by  the 
drag  of  the  mooring  line  itself. 

4.  Assume  that  the  power  at  the  sum  frequency  f  + 
us  is  partially  due  to  the  nonlinear  inter¬ 
action  between  semi-diurnal  tidal  oscillation 
and  a  local  wind-induced  inertial  oscillation. 
Then  shear  at  the  sum  frequency  is  indirectly 
correlated  with  the  wind  at  the  inertial  fre¬ 
quency.  Assume,  in  addition,  that  the  spectral 
shape  of  wind  stress  is  stationary.  Power 
levels  at  all  frequencies  then  rise  and  fall  in 
unison.  Thus  there  will  be  coherence  between 
wind  stress  and  shear  at  the  sum  frequency,  but 
this  coherence  does  not  imply  a  direct  causal 
relationship. 
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Section  5 

CORRELATION  FUNCTIONS 


Auto-correlation  functions,  as  well  as  vertical, 
horizontal,  and  slant  separation  cross-correlation  functions 
were  computed  for  vertical  shear  measurements.  Separate 
calculations  were  performed  for  the  x-  and  y-direction 
components  of  the  shear.  The  measurements  were  not  filtered 
prior  to  computing  the  correlation  functions.  As  a  result, 
spikes  appear  at  zero  lag  for  the  auto-correlations  and  the 
vertical  separation  cross-correlations.  It  was  found  that 
smoothing  the  data  records  with  a  three-point  triangle 
weighting  function  was  sufficient  to  remove  the  spikes.  But 
this  smoothing  was  not  performed  for  the  figures  displayed 
in  the  present  report.  Future  reports  may  examine  the 
appropriateness  of  such  filtering. 

5.1  AUTO-CORRELATIONS  OF  VERTICAL  SHEAR 

Because  the  auto-correlation  functions  of  vertical 
shear  are  all  much  alike,  only  two  example  auto-correlation 
functions  are  displayed  in  Figures  5.1  and  5.2.  Since  auto¬ 
correlation  functions  are  symmetric,  only  the  positive  sides 
are  plotted.  The  spikes  at  zero  lag  are  discussed  above. 
The  strong  ringing  has  a  crest-to-crest  duration  equal  to  the 
inertial  period,  about  19  hours.  The  amplitude  of  the 
ringing  decreases  with  increasing  depth  because  inertial 
energy  becomes  a  smaller  percentage  of  the  total  energy  as 
depth  increases.  This  can  also  be  seen  in  the  frequency 
spectra  shown  in  Figures  3.3  -  3.9.  Such  strongly  periodic 
auto-correlation  functions  near  the  surface  suggest  that  when 


inertial  motions  are  the  dominant  motion,  the  shear  field 
becomes  rather  deterministic.  Thus  a  single  velocity  profile 
measurement  together  with  the  knowledge  of  dominant  inertial 
motions  might  allow  the  reconstruction  at  an  earlier  point  in 
time  of  the  shape  of  an  hydrodynamic  feature  in  order  to 
obtain  clues  as  to  it's  possible  origin. 

5.2  CROSS-CORRELATIONS  OF  HORIZONTALLY,  SLANT,  AND 

VERTICALLY  SEPARATED  VERTICAL  SHEAR 

The  cross-correlations  of  horizontally,  slant, 
and  vertically  separated  vertical  shear  measurements  are 
listed  in  Tables  5.1,  5.2,  and  5.3,  respectively,  and  are 
displayed  in  Figures  5.3  -  5.21.  Since  cross  correlation 
functions  are  not  symmetric,  both  sides  (positive  and 
negative)  sides  are  displayed.  Positive  lags  indicate  that 
the  second  pair  of  vertically  differenced  current  measure¬ 
ments  was  being  lagged  with  respect  to  the  first  pair.  As 
in  the  auto-correlation  functions,  the  ringing  behavior 
appears  to  be  ubiquitous,  with  a  crest-to-crest  duration 
equal  to  the  inertial  period. 

5.3  CROSS-CORRELATION  ANALYSIS 

The  cross-correlation  between  two  sinusoidal  time 
series  is  also  a  sinusoidal  function.  If  the  two  time  series 
are  shifted  in  phase  with  respect  to  one  another,  then  this 
phase  shift  may  readily  be  ascertained  from  the  cross¬ 
correlation  function.  In  this  section  we  present  estimates 
of  the  phase  shifts  shown  by  the  horizontally,  vertically, 
and  slant  separated  cross-correlation  functions. 
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Table  5.1 


CROSS-CORREI.ATIONS  OF  HORIZONTALLY  SEPARATED 
VERTICAL  SHEAR  MEASUREMENTS 


Fig. 

Current  Meters 

Depths  (i 

5.3 

3393-3394/3402-3403 

12-32 

5.4 

3394-3395/3403-3404 

32-52 

5.5 

3395-3396/3404-3406 

52-72 

5.6 

3395-3396/3385-3386 

52-72 

5.7 

3404-3406/3385-3386 

52-72 
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Table  5.2 
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CROSS-CORRELATIONS  OF  SLANT  SEPARATED  VERTICAL 
SHEAR  MEASUREMENTS 


Fig. 

Current  Meters 

Depths  (m) 

5.8 

3393-3394/3403-3404 

12-32/32-52 

5.9 

3393-3394/3404-3406 

12-32/52-72 

5.10 

3394-3395/3402-3403 

32-52/12-32 

5.11 

3394-3395/3404-3406 

32-52/52-72 

5.12 

3395-3396/3402-3403 

52-72/12-32 

5.13 

3395-3396/3403-3404 

52-72/32-52 

5.14 

3385-3386/3393-3394 

52-72/12-32 

5.15 

3385-3386/3394-3395 

52-72/32-52 

5.16 

3385-3386/3402-3403 

52-72/12-32 

5.17 

3385-3386/3403-3404 

52-72/32-52 

% 
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Table  5.3 


CROSS-CORRELATIONS  OF  VERTICALLY  SEPARATED 
VERTICAL  SHEAR  MEASUREMENTS 


Fig. 

Current  Meters 

Depths  (m) 

5.18 

3393-3394/3394-3395 

12-32/32-52 

5.19 

3402-3403/3403-3404 

12-32/32-52 

5.20 

3394-3395/3395-3396 

32-52/52-72 

5.21 

3403-3404/3404-3406 

32-52/52-72 
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(Top)  and  V-Components  (Bottom) . 

(Left):  Between  3395-3396  and  3404-3406. 
(Right):  Between  3395-3396  and  3385-3386. 
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Figure  5.7  Cross-Correlation  Functions  of  Horizontally 
«  Separated  Vertical  Shear, 

U-Component  (Top)  and 
V-Component  (Bottom), 

Between  Current  Meters  3404-3406 
and  3385-3386. 
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Cross-Correlation  Functions  of  Slant  Separated  Vertical  Shear, 
U-Components  (Top)  and 
V-Components  (Bottom). 

Figure  5.16  (Left):  Between  3385-3386  and  3402-3403. 

Figure  5.17  (Right):  Between  3385-3386  and  3403-3404. 
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Cross-Correlation  Functions  of  Vertically  Separated  Vertical  Shear, 
U-Components  (Top)  and 
V-Components  (Bottom). 

Figure  5.20  (Left):  Between  3394-3395  and  3395-3396. 

Figure  5.21  (Right):  Between  3403-3404  and  3404-3406. 
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Figure  5.22  schematically  shows  the  current 
meters,  and  the  cross-correlation  phase  shifts  between  them. 
Only  the  deepest  two  current  meters  from  mooring  338  are 
shown,  as  there  is  no  data  from  the  32  m  current  recorder 
3384  so  there  is  no  shear  time  series  for  the  pair  3383-3384. 
The  lines  join  those  pairs  of  current  meters  for  which  the 
cross-correlation  functions  of  velocity  differ-ences  show 
significant  periodicity  with  an  inertial  period.  Beside  each 
line  is  the  estimated  phase  shift,  in  degrees.  We  estimate 
the  resolution  limit  for  our  phase  shift  estimates  to  be 
about  20°.  The  arrow  through  each  line  points  toward  the 
lagging  pair  of  current  meters.  The  circled  arrows  denote 
significant  asymmetry  in  the  cross-correlation  functions. 

With  but  one  exception  (3394-3395  to  3385-3386), 
the  arrows  for  the  slant  and  vertical  cross-correlations  all 
point  from  the  lower  toward  the  upper  pair  of  current  meters. 
This  result  implies  that  the  phase  velocity  of  these  near- 
inertial  oscillations  is  directed  upward.  From  the  estimated 
phase  shift  we  can  estimate  the  values  of  phase  velocity  and 
vertical  wavelength.  The  phase  shifts  associated  with  the 
vertical  cross-correlations  range  from  60°  to  110°,  cor¬ 
responding  to  upward  phase  velocities  of  1-2  x  10~3 
m/sec.  This  range  corresponds  to  vertical  wavelengths 
between  67  and  134  m. 

The  two  uppermost  horizontal  cross-correlations 
both  show  zero  phase  lag.  Therefore  we  can  only  put  a  lower 
limit  on  the  horizontal  wavelength  in  the  direction  between 
moorings  339  and  340.  Our  estimated  resolution  limit  of 
+20°  implies  that  the  predominant  horizontal  wavelength  is 
at  least  1000  km.  We  remind  the  reader  that  these  cross¬ 
correlation  functions  were  obtained  over  a  50  day  period,  and 
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that  during  portions  of  this  period  the  predominant  hori¬ 
zontal  wavelength  may  be  quite  different. 

Two  of  the  arrows  in  Figure  5.22  are  circled, 
indicating  marked  asymmetry  in  the  correlation  functions 
with  respect  to  2ero  lag.  This  asymmetry  shows  stronger 
inertial  periodicity  for  negative  lags  than  for  positive 
lags.  This  implies  that  the  upper  pair  of  current  meters 
undergoes  near-inertial  oscillations  before  the  lower  pair. 
This  result  agrees  with  our  understanding  that  the  oscill¬ 
ations  are  wind-induced,  and  that  the  direction  of  energy 
propagation  (the  direction  of  the  group  velocity)  is  direct¬ 
ed  downward. 


> 
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SUMMARY  AND  CONCLUSIONS 


We  have  analyzed  vertical  shear  time  series 
obtained  from  48-day  current  meter  records  obtained  by  R. 
Pollard  in  1970.  Three  moorings  in  a  triangular  array  were 
spaced  50  and  70  kilometers  apart.  Each  contained  current 
meters  at  12,  32,  52  and  72  meters  depth.  One  instrument  at 
32  meters  depth  failed  to  give  any  data.  The  moorings  were 
located  on  the  continental  slope  south  of  New  England  at 
depths  of  2400  to  2600  meters.  Two  atmospheric  frontal 
systems  passed  over  the  moorings  near  the  beginning  and  end 
of  the  time  series.  These  fronts  caused  energetic  inertial 
oscillations  which  dominate  the  upper  current  meter  records 
for  several  days  following  the  disturbances. 

Vertical  shear  was  computed  by  simply  taking  the 
difference  between  velocities  recorded  by  vertically  adjacent 
current  meters  and  dividing  by  the  20  meter  sensor  spacing. 
Because  of  the  large  vertical  sensor  spacing,  frequency 
spectra  of  vertical  shear  contain  only  a  portion  of  the  total 
shear  variance.  Based  on  a  spectral  model  of  shear  described 
by  Grabowski  and  Newman  (1981),  we  estimate  that  sensors  20 
meters  apart  in  depth  will  detect  54%  of  the  shear  variance 
which  would  be  measured  by  sensors  10  meters  apart  and  14%  of 
the  shear  variance  which  would  be  measured  by  sensors  spaced 
2  meters  apart.  Sensor  spacing  probably  has  little  effect, 
however,  on  the  observed  shape  of  the  frequency  spectrum 
below  2  cph,  the  highest  frequency  considered  in  this 
analysis.  Therefore,  only  spectral  power  density  levels  are 
affected. 


Frequency  spectra  of  all  the  shear  time  series 
contain  strong  peaks  in  the  neighborhood  of  the  inertial 
frequency  (0.053  cph  or  9.3  x  10” 5  rad/sec).  The  variance 
associated  with  the  inertial  peaks  decreases  with  increasing 
depth.  Vertical  shear  at  inertial  frequencies  is  created 
largely  by  the  rapid  decay  of  velocity  amplitude  with 
distance  from  the  surface.  That  is,  inertial  velocities 
recorded  at  the  upper  end  of  a  given  20  meter  interval  are 
substantially  greater  that  those  recorded  at  the  lower  end. 

Frequency  spectra  of  shear  also  show  a  distinct 
peak  at  the  semidiurnal  frequency,  but  such  peaks  typically 
contain  substantially  less  energy  than  an  inertial  peak  in 
the  same  spectrum.  The  peak  near  the  semidiurnal  frequency 
is  present  because  the  semidiurnal  tide  is  strongly  baro- 
clinic  in  the  New  England  slope  region.  According  to  Regal 
and  Wunsch  (1973),  the  baroclinic  semidiurnal  tide  is  due  to 
the  interaction  of  the  barotropic  tide  with  the  continental 
slope.  However,  the  semidiurnal  tide  is  not  expected  to 
contribute  significantly  to  frequency  spectra  of  near  surface 
shear  in  the  open  ocean. 

Shear  spectral  density  levels  were  found  to 
increase  with  frequency  in  a  band  from  0.5  to  2.0  cph.  This 
effect  may  be  due  to  aliasing  of  the  internal  wave  motions  at 
frequencies  higher  than  2.0  cph.  An  example  of  the  aliasing 
effect  was  discussed  based  on  the  Garrett  and  Munk  internal 
wave  spectrum.  Mooring  motion  may  be  another  factor  contri¬ 
buting  to  the  spectral  level  increase  with  frequency. 


Evidence  of  a  possible  nonlinear  interaction 
between  motions  at  inertial  and  semidiurnal  frequencies  was 
found  in  a  spectral  peak  at  a  frequency  equal  to  the  sun  of 
the  inertial  and  semidiurnal  frequencies.  Wave  motions  at 
this  sum  frequency  can  propagate  into  the  stratified  interior 
and  thereby  provide  radiative  damping  of  inertial  motions. 
Other  suggested  dissipation  mechanisms  are  Kelvin-Helmholtz 
instability  (Pollard  et  al. ,  1973)  and  radiating  internal 
waves  generated  by  differentially  advected  mixed  layer 
turbulence  inducing  undulation  of  the  seasonal  thermocline 
(Bell,  1978). 

Because  of  the  failure  of  one  32m  current  meter, 
horizontal  coherence  was  computed  for  only  one  pair  of  shear 
time  series  in  the  12-32  meter  range.  This  horizontal 
coherence  spectrum  (over  a  distance  of  50  kilometers)  showed 
significant  coherence  in  the  frequency  range,  -f<u><0.  The 
corresponding  phase  spectrum  fluctuates  about  zero  in  the 
same  frequency  range  implying  a  horizontal  wavelength  very 
long  compared  to  50  kilometers. 

Coherence  of  shear  with  wind  stress  was  significant 
at  approximately  0.7  in  the  band  -f<u><0,  with  an  expected 
resonant  phase  shift  of  180°  near  u>  *  -f.  In  the  depth 
intervals  32-52  meters  and  52-72  meters,  almost  no  signifi¬ 
cant  coherence  was  found  between  shear  and  wind  stress.  The 
latter  observation  is  consistent  with  the  propagation  of 
near-inertial  oscillations  into  these  depth  intervals  along 
very  gently  downward  sloping  characteristics  from  a  remote 
surface  location,  forced  by  the  wind  stress  at  an  earlier 
time.  Some  similarities  were  observed  between  these  results 
and  those  of  Weller  (1981)  derived  from  data  taken  by 
profiling  current  meters  from  FLIP. 
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Some  peaks  in  conerence  between  wind  stress  and 
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2)  effects  of  imperfect  parameterization  in  calculating  wind 
stress,  or  3)  effects  due  to  correlation  between  changing 
power  levels  of  a  constant  shape  wind  spectrum  and  the 
corresponding  changes  in  level  of  inertial-tidal  inter¬ 
action. 

Vertical,  horizontal,  and  slant  cross-correlation 
functions  were  computed  for  the  vertical  shear  time  series. 
All  of  the  vertical  cross  correlation  functions  and  all  but 
one  of  the  slant  cross  correlation  functions  show  upward 
phase  propagation  corresponding  to  phase  velocities  in  the 
range  1-2  x  10-3  m/sec,  and  vertical  wavelengths  in  the 
range  67-134  meters.  Two  of  the  cross-correlation  functions 
have  amplitudes  that  are  larger  for  negative  than  for  posi¬ 
tive  lags.  The  implication  is  that  internal  wave  packets 
propagate  downward,  away  from  the  near-surface  region  of 
generation. 

Several  conclusions  may  be  drawn  regarding  the  open 
ocean  near  surface  shear  observed  in  this  data  set.  First, 
because  the  data  were  obtained  in  deep  water  well  away  from 
land,  they  should  represent  open  ocean  conditions  rather  well 
with  the  one  exception  of  the  shear  caused  by  the  strongly 
baroclinic  semidiurnal  tide.  In  the  open  ocean  the  semi¬ 
diurnal  tide  is  barotropic  and  cannot  produce  shear.  Second, 
it  is  clear  from  these  data  that  inertial  motions  can  cause  a 
large  fraction  of  the  shear  calculated  over  a  vertical  scale 
of  20  meters.  Consequently,  the  vertical  shear  field 


exhibits  some  horizontal  and  vertical  scales  associated  with 
inertial  motions.  Third,  since  dissipation  tines  for 


iner 

tial  motions 

are  of  the  order 

of  several  days, 

and 

typical  intervals 

between  passages 

of  atmospheric  f 

rents 

may 

be  comparable 

in  northern  areas 

,  then  inertial  mo 

tions 

may 

exist  much  of 

the  time.  This 

implies  that  the 

shear 

field  may  have  a  large  inertial  component  and  therefore  in 
fact  be  quasi-periodic  much  of  the  time.  Fourth,  fifteen 
minute  sample  intervals  are  too  long  and  cause  aliasing  of 
the  internal  wave  induced  shear.  The  shorter  112.5  second 
sample  intervals  used  on  Vector  Averaging  Current  Meters 
( VACMs )  are  much  more  likely  to  sample  all  important  internal 
wave  frequencies  adequately.  Finally,  although  some  effects 
of  mooring  motion  may  be  subtracted  out  when  differencing 
velocity  time  series,  coherence  peaks  at  w  *  -(f  +  u^)  are  an 
indication  of  possible  effects  which  remain.  Analyses  of 
shear  obtained  from  current  meters  must  therefore  be  done 
with  some  caution  to  avoid  erroneous  results  due  to  these 
effects. 
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